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Abstract 

The invariant mass spectrum of the e~^e~ and f^~^fJ-~ pairs from decays of 
light unflavored mesons with masses below the (/)(1020)-meson mass to final 
states containing along with a dilepton pair one photon, one meson, and two 
mesons are calculated within the framework of the effective meson theory. 
The results can be used for simulations of the dilepton spectra in heavy-ion 
collisions and for experimental searches of dilepton meson decays. 
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I. INTRODUCTION 



During the last decade, the problem of the description of hadrons in dense and hot nuclear 
matter received great attention. It is well known that particles change their properties when 
they are placed into a medium. Already in the seventies the reduction of the nucleon masses 
in nuclei was implemented into the Walecka model in the framework of effective hadron 
field theory Later on this effect was put on firmer grounds on the basis of a partial 

restoration of chiral symmetry and finite-density QCD sum rules The change of the 
meson properties is also discussed in the quantum hadrodynamics 0] and finite-density QCD 
sum rules [Q. The investigations in the Nambu-Jona-Lasinio model provide an evidence for 
reduction of nucleon and meson masses at finite density and temperature as well |^. 

On the other hand, many-body correlations lead to a dressing of the particles inside the 
medium and strongly modify spectral properties of the mesons Thus one expects a 

significant reduction of the corresponding life times which, loosely speaking, result in melting 
the mesons in the nuclear environment. 

The melting of the higher nucleon resonances is established experimentally from the 
measurement of the total photoabsorption cross section on heavy nuclei . There is a clear 
signal for a change of the shape of the A-resonance in nuclei, related to the Fermi motion 
and, as noticed in Ref . , to the collision broadening effect discussed first by Weisskopf |TT 



in connection with a broadening of the atomic spectral lines in gases. The higher nucleon 
resonances are not seen in the cross section due to a strong collision broadening effect. 

The purpose of the current investigations is to determine mass shifts and the broadening 
of hadronic resonances in nuclear matter. The formulation of this problem can be traced 
back to the atomic spectroscopy where the shifts of atomic energy levels and the broadening 
of atomic spectral lines in dense and hot gases is a relatively well studied subject (see e.g. 
Ref. PI). 

The search for signatures of modified hadron properties such as reduction masses ac- 



cording to the Brown-Rho scaling |T3[ are presently pursued with high experimental efforts. 
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However, most hadronic probes loose important information on the early and violent phase 
of the reaction by strong final state interactions. One of the most promising probes to 
study in-medium properties of vector mesons are dilepton pairs [i = e, fi) which are 
produced from decays of p^-, uj-, and 0-mesons in heavy-ion collisions. The leptonic probes 
have the advantage that they are the nearly undistorted messengers from the conditions at 
their creation, unlike e.g. pions. 

The dilepton spectra measured by the CERES and HELIOS-3 Collaborations at CERN 
SPS have attracted special interest 0,^. Compared to theoretical predictions, both ex- 
periments found a significant enhancement of the low-energy dilepton yield below the p 
and u peaks. One way to explain this low-energy dilepton excess is to assume the scenario 



of a significant reduction of the p- meson mass in a dense medium [jT6|-|T^. On the other 
hand, a more sophisticated treatment of the p-meson spectral function which includes the 
broadening in dense matter seems also to be sufficient to account for these data [|18|,|19 



An excess of low-energy dileptons occurs already at moderate bombarding energies. How- 
ever, the spectra obtained by the DLS Collaboration at the BEVALAC for the incident 



energies around 1 A-GeV cannot be reproduced by present transport calculations |2T|]. Even 
using the reduction of the p-meson mass, there remains a discrepancy for the dilepton yield 
by a factor of 2 to 3 ||2l|j2^. Also the medium dependence of the spectral function, even 
in combination with a dropping p mass, does not provide an explanation for this so called 
'DLS puzzle' [p3| , p4| . It is interesting to note that this fact is independent on the system size 



and occurs in light {d + Ca) as well as in heavier systems {Ca + Ca). The dilepton yield in 
elementary p + p collisions measured by the DLS Collaboration p5[ is also underestimated 



by standard theoretical descriptions In Ref. it was claimed that the discrepancy 
between the DLS data and theoretical simulations disappears if additional background con- 
tributions from dilepton decays of higher nucleon resonances (mainly A^*(1520)) are taken 
into account. This demonstrates that a most precise knowledge of the background is in- 
dispensable for the interpretation of the present and future dilepton data. The HADES 
experiment at GSI, Germany, will focus on these topics to a large extent |]27|| . 



While the dilepton spectra stemming from the mesonic decays are not distorted by final 
state interactions, the problem of extracting information on the in-medium properties of the 
vector mesons is a specific theoretical task, due to a large amount of decays contributing 
to the background. A precise and rather complete knowledge of the relative weights for 
existing decay channels is therefore indispensable in order to draw reliable conclusions from 
dilepton spectra. 

The study of the dilepton decays is useful also for the search on the dilepton decay modes 
of the light unfiavored mesons. These decays give a deeper insight into meson structure, al- 
lowing to measure transition form factors at the time-like (resonance) region. In four-body 
decays the decay probabilities are determined by the half-off-shell meson transition form 
factors which cannot be measured in other reactions. There are plans to study the r^' tran- 
sition form factors in the space-like region in reactions of the photo- and electroproduction 
at CEBAF energies [|2^. Recently, the DLS Collaboration published data on the dilepton 



production in the elementary pp and pd collisions [^. These data are analyzed in Refs. 
PT| , |2B[] . The plans from the HADES Collaboration include measurements of the dilepton 
spectra from proton-proton and pion-proton collisions ET^. These experiments are stimu- 



lated by the already mentioned discrepancy between the number of the observed dilepton 
events and results of the transport simulations, that indicate a limited understanding of the 
mechanism for the dilepton emission in heavy-ion collisions. In such experiments, there is a 
possibility for exclusive measurements of the different mesonic channels. 

In the present work, we perform a detailed study of possible mesonic decays 
which appear in the SIS energy range, i.e. at lab. energies below 2 GeV/nucleon. At 2 
GeV/nucleon, the p-meson is produced slightly above the threshold. Due to statistical fluc- 
tuations, the production of heavier mesons is, however, also possible. We consider decays of 
unflavored light mesons with masses below the 0(lO2O)-meson mass within the framework of 
the effective meson theory. The vertex couplings are determined from the measured strong 
and radiative decay widths and, when the experimental data are not available, from SU{3) 
symmetry. The transition form factors entering the decay rates are calculated using the 



Vector Meson Dominance (VMD) model. In this way, we achieve generally good agreement 
with the experimental branching ratios for radiative meson decays. For the dilepton decay 
modes the branching ratios are, however, only known in a few cases. 

We consider the vector mesons p, u, and 0(1020) (= V), the pseudo-scalar mesons tt, 
7^, r^' (= P), and the scalar mesons /o(980) and ao(980) (= S). The various decay modes 
can systematically be classified as follows: (i) There are, first of all, the direct decays modes 
V — > £'^£~, which contain the information on the in-medium vector meson masses. In the 
next Sect., some useful relations which simplify calculations of the decay rates to final states 
with a dilepton pair are derived. In Sect. 3 we make a few remarks on the direct decay modes. 
There exits then a large number of processes which mask the vector meson peaks and which 
should be treated as a background, (ii) These are Dalitz decays of pseudoscalar mesons 
P — >• ■yi'^i~ and scalar mesons S —>■ . These decays are discussed in Sect. 4. (iii) One 

has also Dalitz decays with one meson in the final states V — > Pi'^i~ and P — > Vi'^i", 
which we discuss in Sect. 5. The radiative decays of these mesons are well studied both 
from the theoretical and experimental points of view. The uncertainties in the estimates 
for the dilepton decays are connected with the lack of experimental information on the 
transition form factors. Constructing the VMD model transition form factors, we take 
special care of the quark counting rules, (iv) Finally, there exist decays to four-body final 
states V PPe+£-, P PP£+£-, and S PPe+l'. Sect. 6 is devoted to these decays. 
We calculate almost all four-body dilepton modes of the unflavored mesons with masses 
below the 0(1020). The numerical results for radiative widths of the three-body decays 
{V PPj, etc.), which provide a useful test for the model considered, for the dilepton 
widths {V — > Pi'^i~, V — > PPt^l~, etc.), and for the dilepton spectra from the unflavored 
meson decays are presented in Sect. 7. 
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II. RELATION BETWEEN THE DECAYS M ^ M'7* AND M ^ M'l+i' 



As mentioned above, we consider decays M M'i~^i~ where M is a meson, M' is a 
photon, a meson, or two mesons, and £~^£~ is an electron-positron or muon-antimuon pair. 
The results of this Sect, are valid, however, for arbitrary states M'. The decay M M'i^i~ 
proceeds through two steps: M M'7* and 7* , where 7* is a virtual photon whose 

mass, M, is equal to the invariant mass of the dilepton pair. 

The matrix element of the physical process M M'7 for a real photon 7 has the form 

M = M,e;{k) (II.l) 

where s^lk) is a photon polarization vector. The matrix element Ai fj, is defined also at 
k"^ = 7^ for virtual photons 7*. As a consequence of the gauge invariance, it is 
transverse with respect to the photon momentum 

Mf,k^ = 0. (II.2) 

The decay rate M M'7* can formally be calculated as 

dm M'7*) = ^E^H>v(-^7..)^^j^rf^^^^ (II.3) 

where ^/s is mass of the decaying meson and n is number of particles in the state M'. The 
phase space in Eq. (|II.3|) is defined in the usual way 

d^,{^, m„ m,) = n ^,S\P - E (11-4) 

i=l * i=l 

Here, P is the four-momentum of the meson M , = s, and Pi are momenta of the particles 
in the final state, including the virtual photon 7*. In Eq. (|II.3|) , the summation over the final 
states and averaging over the initial states of the decaying meson is performed. The limit 
gives the decay rate of the physical process M — > M'7. 
The M — > M'i^i~ decay rate is given by 

dm ^ = ^ E^H>^7j.j/^^|^ (II.5) 
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where is the lepton current. The term comes from the photon propagator, and 

is the phase space of n particles in the state M' and of the pair. 
The value r(M M'i+i-) can be related to the decay rates r(M M'7*) and 1(7* 
In the analogy with massive vector particles, the width of a virtual photon 7* can 
formally be evaluated as 



Mr(7* - tn = -{M' + 2mj)^l - (II.6) 

where is the lepton mass. The expression for the product of two dilepton currents, 
summed up over the final states of the i'^i" pair, has the form 

Y.^,J.* = ^(M^ + 2mj)i-g,. + ^) (II.7) 
/ 

where a is the fine-structure constant and k is the total momentum of the pair. Factorizing 
the ra-body invariant phase space, 

rf$fc(v^, mi, ruk) = rf$A,._i(v^, mi, mfc_2, M)dM'^^2{M, mfc_i, m^), (II.8) 

which can be proved by inserting the unity decomposition 

1 = J d^qdMH{q^ - M^)5\q - pu-i - Pk) (II.9) 

into Eq.(|L|), one obtains from Eqs.(|L3]), and (|L6D with the help of Eqs. (p7D and 

( p.8| ) the following expression 

dT{U M'£+r) = dT{U M'7*)Mr(7* ^ ^^^")^^- (^^-^O) 

The factor dM"^ j (vrM^) has the form of a properly normalized Breit-Wigner distribution for 
a zero-mass resonance. 

The two-body phase space in Eq. ( [II. 8|) has the form 

<l>2(Vs,mi,m2) = p (11.11) 

where 



p (Vs,mi,m2) = ^-^= (11.12) 

is momentum of the particles 1 and 2 in the c. m. frame. 

In the following, we will work with the matrix elements of the processes M — > M'7* and 
use Eg. ( [II. 3D to derive results for the decay rates M — > M'7 with real photons and Eq.( [II.10D 
to get results for decay rates M V^' t^i' with dileptons in the final states. Now, we will 
consider the processes which are interesting for the study of dilepton spectra in heavy-ion 
collisions. 



III. DECAYS OF THE p-, w-, AND 0-MESONS TO PAIRS 

The diagram for the V decays with ^ = p, cj, and (j) is shown in Fig. In terms 



of the vector meson fields, V^, the electromagnetic current has the form [29 



3, = -eY.—y^ (III-l) 

9v 



A" 

V 

where my are the vector meson masses and e = — |e| is the electron charge. The SU{'i) 
predictions for the coupling constants, gp : gu, : g^j) = I : 3 : are in good agreement with 
the ratios between the values gp = 5.03, g^ = 17.1, and g^p = —12.9 extracted from the e^e~ 
decay widths of the p-, u-, and 0-mesons with the use of the well known expression 

Yiy ^ tl-) = -^(1 + 2^y(mv, m,, m,), (III.2) 

3gv '^v 

with p* defined in Eq. (|II.12|) . 



IV. MESON DECAYS TO PHOTONS AND PAIRS 

In this Sect., we discuss meson decays P ^t^t~ (Fig. 2) and 5* — > ^t^i" (Fig. 3) 
for P = 7r°,?7, and rj' and S = /o(980) and ao(980). As we shall see, these decays are the 
dominant e^e^ modes for vr*^-, ?7-mesons, for r^'-meson at M > 250 MeV, /o-meson at M 
> 500 MeV and in some other cases. The fi~^fi~ modes are also discussed. 
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These decays are related to the experimentally measured two photon decays. The un- 
certainties in the estimates originate only from the purely known transition form factors in 
the time-like region. The ^777* transition form factor is in reasonable agreement with the 
one-pole VMD model predictions The experimental errors in the r( transition form 



factor are large The one-pole VMD approximation for the P77* transition form factors 
is in agreement with the quark counting rules which predict for these form factors a ~ 1/t 



asymptotics |32 



The nature of the scalar mesons has been a subject of intensive discussions for a long 
time. The 4-quark content of the scalar mesons would imply a ~ \jt^ asymptotics for the 
5*77* transition form factors. In order to provide the correct asymptotics for the 4-quark 
meson transition form factors, the VMD model should be extended to include contributions 
from higher vector meson resonances. The 07717 transition form factor has also ~ 
asymptotic behavior |3^. This form factor is measured in the time-like region and the 



data show deviations from the naive one-pole approximation. The inclusion of higher vector 
meson resonances improves the agreement and provides the correct asymptotics. 

The results of a recent measurement of the 7/0 branching ratio are interpreted as 
evidence for the dominance of a 4-quark MIT bag component in the /o-meson wave function 



3J,|35|. We thus calculate branching ratios S 'yi'^i" assuming the 4-quark nature of the 
scalar mesons and imposing constraints from the quark counting rules to the 6*77* transition 
form factors. 

A. Decay modes vr*^ — > 7e"'"e~, r] 'y£~^£~, and 77' — > 'yi~^£~ 

The effective vertex for the P — 77 decays has the form 

dMuA^P (IV.l) 

where P = n^, r], and t]' and is the photon field. The matrix element for the decay P 
77* with a virtual photon 7* has the form 
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M = -ifp^^Fp^^{M'')er.^,Kel{k)h^el{h) (IV.2) 

where k is the virtual photon momentum (/c^ = M^), ki is the real photon momentum 
{kl = 0), and Fp^^{t) is transition form factor P77*. The comparison of the P — > ^t^t~ 



decay width with the decay width of a physical process P 77 allows to write [30 



2 K \^ ' ' M Fp^ (M') Mr(7* -^£+r)— — (IV.3) 

where i/i = /xp is the pseudoscalar meson mass. The value Mr(7* — s> is given by 

Eg. ( [11.61) . The factor 2 in Eq.( |lV.3[) occurs due the identity of photons in the decay P — 77. 
The product of the qubic term and the absolute square of the form factor gives the ratio 
between squares of the matrix element (P^V.2|) at k^ = and k"^ = 0, multiplied by the ratio 



between the two-particle phase spaces. The quark counting rules [^] imply that the form 
factor Fp^^(t) behaves as~l/)f:att^oo. The experimental data are described reasonably 
well by the monopole formula 

Fp,,{t) = (IV.4) 

with Ap = 0.75 ± 0.03, 0.77 ± 0.04, and 0.81 ± 0.04, respectively, for the 7r°-, r]-, and r]'- 
mesons which reproduces the correct asymptotics. Such a monopole fit can naturally 
be interpreted in terms of the vector meson dominance. The values of the Ap's are close 
to the p- and cj-meson masses. In a QCD interpolation formula, Ref. |3^, the pole masses 



are related to the PCAC coupling constants of the pseudoscalar mesons, Ap = 27r/p. This 
expression gives for the pole masses similar numbers. 

In the T]' — i> ■ye~^e~ and t]' — >• decays, the p° and u poles occur in the physical 

region allowed for the spectrum of the dilepton pairs. The SU (3) symmetry with the 77 — 77' 
mixing angle 9 = arcsin(— |) = —19.5° which is quite close to the experimental value 
Qexp _ _i5 5 -|- 13 (^ggg g_(^_ i^g) predicts /p^^' = Sf^^r}' = :^/a;77r- The experimental 
branching ratios P"^P(r/' 7p°) = (30.2 ± 1.3)% and P"^P(r/' ju) = (3.02 ± 0.30)% 
indicate that the ratio between the p°- and cu-meson couplings with the r/'-meson equals 
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fp-rri' / fu-yr}' — good agreement with the SU{3) predictions. In the timehke resonance 

region, the transition form factor ( [IV .41 ) should be modified: 

Fn',iit) = Cp , ' . ^ + , ^ + cx^^ (IV.5) 

— t — tnipL p — t — tm^jL uj f^x — t 

(recall that A^/ ^ nip ^ m^), with the weight coefficients Cp ~ fp-yri'/gp = fppri'/gp and 
Cu ~ fivw/gu, = fuujTi'/gl- The SU{3) symmetry predicts fppr,> = f^^r,' = ^fp^-n- The 
normalization condition looks like Cp+ c^; = 1. The slope is close to From the 

relation gp/guj ~ 0. 3 one gets Cp/cuj ~ 10.5, the relative sign of the values Cp and c^; is 
fixed by SU (3) symmetry. The third term in Eq. ([IV. 5|) is introduced for reasons explained 



below. We assume for the moment cx = 0. With Cp = 0.9 and Cuj = 0.1, we obtain 
B{ri' — s> 7yU+yU^) = 0.90 x 10"^, in good agreement with the measured value B'^^^ij]' — >■ 
= (1.04 ± 0.26) X 10^^. Data for the rj' e"^e~7 decay are not available. The 
coupling constants we used in the calculations are summarized in Table 1. 

The product of the branching ratios B^^^{r]' '~fp^)B^^^{p^ —>■ p,^ = 1.4 x 10~^ is 
almost one order of magnitude smaller than the value B^^'^ij]' —>■ The magnitude of 

the direct p°-meson contribution to the rj' '~ffi~^fi~ decay can be estimated from Eq.( [lV.3] ) 
with the use of a narrow-width approximation for the transition form factor, 

\F,,,,{t)\'^\cpf—^S{t-ml), (IV.6) 
p 

to give B{7]' 7p° — > 7yU+yU^) ^ 2 x 10~^. This value is already close to the product 
B^^P[r]' —>■ p^) B^^^ {p^ p^p^), but still higher. The relative contributions of the p"- 
and tu-mesons to the 77' — 'ye~^e~ and rj' 'yp~^p~ decay rates are inversely proportional 
to the vector meson widths, as it follows from Eq. ( [IV. 6]) . Since the tu-meson width is only 
8.5 MeV, its contribution is strongly enhanced. The direct contribution of the cu-meson 
to the transition t]' —>■ 'yp'^p~ equals B^^^{ri' —>■ '~fuj)B{uj —>■ p^ p~) ~ 2 x 10~^ where 
B{uj —>■ p^p^) ~ B'^^^{uj e^e^) = (7.15 ± 0.19) x 10~^ (the quoted experimental values 
are all from Ref. |^). The use of the narrow- width approximation gives B{ri' —* ^uo — >■ 



p'^p 7) ~ Fp/r^ \cuj/cpf B{r]' — > 7p° — i> 7p^p ) 4 x 10 ^. This value is also greater than 
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the product B^^^irj' — > ^uj)B{uj — > It can be interpreted as a noticeable contribution 

from the excited vector mesons to the form factor 

In any case, the above estimates demonstrate that the direct and a;-meson contribu- 
tions to the 77' — > 7//+//" decay are small. Thus, the value B{r]' — > 7//"^//") is determined 
mainly by the background, i.e. by those values of which are not near to the p°- and 
a;-meson poles. Since we are interested in the spectrum of the dilepton pairs, it is impor- 
tant, however, to fix the relative weights of the vector meson contributions. A 10 — 20% 
decrease of the residues Cp and c^j due to the admixture of an excited vector meson X yields 
for the measured decays 77' — > j/i'^iJ,', rj' — > 7p°, and 77' — > 70; a consistent interpretation. 
With Cp = 0.8, = 0.08, and cx = 0.12 we get ^(77' -fp° 7A«+A«~) ~ 1-5 x 10"^ and 
B(rj' ^ijj ^ Ri 2 X 10~^, in good agreement with the products of the branching 

ratios B'^^'^{r]' ^p^)B^'^'^{p^ A*"^/^") and B'^'^'^{r]' ^uj)B{uj — >■ p'^p~). The mass of 
the radially excited vector meson X is not well fixed (see discussion in Sect. 5). We set 
mx — 1.2 GeV. Since it is out of the physical region, the width Tx is set equal to zero. 

B. Decay modes /o(980) ^ 7^+^" and ag(980) ^ 7^+^" 

The isoscalar /o(980)-meson and the isotriplet ao(980)-meson have quantum numbers 
jG^jPC^ _ o+(0++) and 1~(0+'^). The /o- and the neutral oo-mesons decay to two photons. 

The effective vertex for the 5* — > 77 decay has the form 

6Cs^^ = fs-i-iFr^FT-pS (IV. 7) 

where — dp At — drAp. The matrix element for the process S — > 77* is given by 

M = -ifsjjFs^'yiM^){9Tahx - 9Txha)ignakx " gpxka)£*riki)e*p{k) . (IV.8) 

Here, ki and k are real {kf — 0) and virtual (/c^ = M^) photon momenta, Fs^^{t) is transition 
form factor of the decay S 77*. The square of the matrix element summed up over the 
photon polarizations equals 8sp*^{^/s,0, M), with ^/s = nis being the scalar meson mass. 
The width of the 5" je'^e~ decay can be written as follows 
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ns - 77) " P*'{v~s,o,o) 1^"^^^^ ^1 ^^^^ ^ ^ ^ W- ^^^-^^ 

The value of r(/o — > 77) is given by PDG [^. In case of the ao-nieson, PDG gives the 
quantity r(ao 77)r(ao TT'^ri)/Ttot{o,o)- Since the mode ao vr^r/ is the dominant one, 
the two-photon width can also be estimated. 

The transition form factor Fs^yyit) depends on the nature of the scalar meson S. The 



asymptotics of the form factor according to the quark counting rules is ~ 1/t for a 



2-quark model and ~ for a 4-quark MIT bag and a KK molecular models of 



the /o" and ao-mesons |^2|. As pointed out by N. Achasov and Ivanchenko the decay 
TT^Ti^'-^ can provide an important information on the structure of the /o-meson. In 
the SND experiment at the VEPP-2M e~^e~ collider the branching ratio 7r°7r°7 was 



measured p^ , |35| . The decay goes mainly through the 7/0 decay mode. The results 
are consistent with the hypothesis of a 4-quark MIT bag nature of the /o-meson. We thus 
calculate dilepton branching ratios S ^£'^i~ assuming the 4-quark nature of the scalar 
mesons. 

In the framework of the VMD model, the 6*77* transition form factor can be reproduced 
assuming contributions from ground-state and excited vector mesons with masses rrii 

i * 

The normalization F5^^(0) = 1 and the asymptotic condition Fs.y.y{t) ~ at t cxd give 
constraints to the residues q: 

l = J2c,, (IV.ll) 

i 

= ^Ciml (IV.12) 

i 

In case of the form factor /o77*, the OZI rule implies that uu, pp and 00 contributions 
are small, while the contributions up, 0p are forbidden by isospin conservation. The relative 
weights of the uj- and -mesons in the form factor are fixed and thus at least three vector 
mesons should be considered to fit the asymptotic behavior: 
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The overall normalization factor can be derived from the condition Fy„^^(0) = 1, the value 
cx is fixed from the requirement Ff^^^it) ~ as i — > oo. The third meson mass mx is 
not well determined. The resulting form factor is given by 

m'^m'^mx{l + Ct) 
= {mi-t){ml-t){mj,-t) ^^^"^^^ 

where 

mliml ~ ml) + mlimj^ - ml) i ^^ 

^ ~ 2 27^ 2 2 2l ■ (^iV.iOj 

The transition form factor a^^yj* has the same structure with the replacement cu p. 

The virtual photon from the scalar meson decay lies in the physical region of the oo- and 
p-mesons, so the cu- and p-meson propagators should be modified by introducing the finite 
meson widths. We use mx — 1.2 GeV. The pole t — mx is outside of the physical region, 
so we set Fx — 0. 

V. MESON DECAYS TO ONE MESON AND A PAIR 

The radiative decays V — > P7 and P — > are well studied experimentally. The 
dilepton modes are measured only for decays 00 n^i^i^ and r]e~^e~ . We use data on 
the radiative decays to fix coupling constants entering the VPj vertexes. These values are in 
good agreement with the SU{3) symmetry relations. We use the framework of the extended 
VMD model and impose constraints to the residues of the transition form factors from the 
quark counting rules. The calculation of the dilepton spectra is then a straightforward task. 

A. Decay modes u — Tr^e+e", p — Tre+e", and cf) — 'K^e^e~ 

The effective vertex for the V P7 radiative decay with V = a;, p°, and (p and P = 7r°, 
77, and rj' has the form 
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SCv^p = -efv-yperafiudaVrduA^P. (V.l) 

The matrix element for the V — > P7* transition can be written as 

M = -iefv^pFv^p{M^)e^^^,tr{P)P.el{k)K. (V.2) 

The diagram for this transition is shown in Fig. ^ The vertex form factor Fy-ypit) normal- 
ized at t = to unity depends on the square of the photon four-momentum, t = M^. The 
width of the decay V — > P7* has the form (see. e.g. p4| , ^ ) 



r(y ^ PY) = \Fv,p{k')\'p*%V~s,ftP,M) (V.3) 

where y/s = my is the vector meson mass and /ip is the pseudoscalar meson mass. The 
limit = describes the radiative decay V P7. Using the experimental value of the 
oj — > 7r°7 width, one gets fui-y-w = 2.3 GeV""*^. The decay width V — > Pi^i^ is connected to 
the radiative width V P7 [Q: 

where Mr(7* t^l~^ is given by Eq.( [11.6|) . The first factor is the ratio between squares of 
the matrix element ([V.2| ) at A;^ = and A;^ = 0, multiplied by the ratio of the two-particle 
phase spaces $2(v^, /^p, ) and $2(v^, /Wp, 0). 

The vertex form factor Fy.yp{t) falls off asymptotically as I^Bj. The additional power 
1/t as compared to the asymptotics of the pion form factor occurs because of suppression 
of the quark spin-flip amplitude due to conservation of the quark helicity in interactions of 
quarks with gluons. 

The wyiT transition form factor can be reproduced taking the contribution of the ground- 
state p- meson and at least one excited p-meson into account. The extended VMD model 
with rrip = 769 MeV and mp> = 1450 MeV provides the correct asymptotics and yields a 
better description of the experimental data p7| , |48[| than the one-pole model with only the 
ground-state p-meson. However, it still underestimates the slope of the form factor at t = 0. 
The experimental value P;^^^(0) = 2.4 ±0.2 GeV"^ prefers a value mp> ^ 1.2 GeV. 
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The possible existence of vector mesons / = 0, 1 with masses around 1.2 GeV was 
under discussion for a long time. A quite strong evidence for significant contributions to the 
spectral functions of the nucleon form factors at values t lower than 1.45 GeV comes from 
the fact that the experimental data for the Sachs form factors are reasonably described by 
the dipole formula 

GEp{t) ^ GMp{t)/fip ~ GMn{t)/f^n ~ _ (V.5) 

where t is in units GeV^. This formula has a double pole. The spectral functions are 
proportional to 

6{t - ml) - 6{t - m\) 



(5'(t-0.71^ 



The double pole in Eq.(V.5) indicates an enhancement in the spectral functions at mx close 
to mp, whose nature is not yet clear as long as the existence of the non-strange vector mesons 
at 1.2 GeV is not established. In any case, it is desirable to shift masses of the first excited 
non-strange "vector mesons" to lower values. This has been done in some popular fits of 
the nucleon form factors (see e.g. [0). The more recent fits use, however, higher values of 
the vector meson masses |Q. In our calculations, the value mx is assumed to be 1.2 GeV, 
as extracted form the slope of the F^^T,{t) at t = 0. 

The multiplicative representation of the transition form factor 07771, which in the zero- 
width limit is completely equivalent to the additive representation, is given by 

2 2 
m,m-Y 

The experimental data on the transition form factors are available for the ?7'77*, ''777*, 
and ti;7r7* transitions. The parametrization (IV. 5) for the Vj form factor is in good agreement 
with the data from Ref. [^, Fig. 26. For the //-meson form factor (IV. 4), our curve coincides 
with the VMD curve from Fig. 24, Ref. [^, since the parameter A = 0.77 GeV in Eq.(IV.4) 
is close to the p-meson mass. 

The one-pole VMD model is known to be in rather pure agreement with the data on the 



iij-K^* transition form factor The two-pole extended VMD model, Eq.(V.6), improves 
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the agreement without introducing new parameters. It fits well the experimental points at 
t < 0.2 GeV^, but underestimates the last three ones from Ref. below 0.4 GeV^. By 
price of introducing a new vector meson, p", and one fitting parameter, the last three points 
can be described better. The modification consists in inserting an additional multiplier 
m^„(l + Ct) I ivn}p„ — t) into the right side of Eq.(V.6), with C being a free parameter. This 
model again satisfies the quark counting rules. However, it overestimates the form factors 
values at t ~ 0.2 GeV^. For the present calculations, we stay in the framework of the 
simplest two-pole VMD model, Eq.(V.6). It is clear, however, that if one needs to calculate 
the dilepton production for the DLS and/or future HADES data, the sensitivity of the results 
on the UTT'y* form factor values must be investigated. 

The 0- meson decay width is given by Eq. ( [V.4| ). In the 0-meson decay, the emitted photon 



has isospin J = 1, so the transition form factor has the same form as for the cj-meson. The 
i~^£~ pairs appear in the physical region of the decay 7r°p°. In order to include the 
p-meson direct contribution, the form factor ( |V.6| ) should be modified by introducing the 
finite p-meson width. In the cases of the u vr^e^e^and p — » vre+e^ decays, vector mesons 
appear in the unphysical regions and their widths are not so important. 

The direct contribution to the 7r°p° decay can be estimated with the use of the 
narrow-width approximation for the transition form factor in analogy with the decays t]' — > 
'-fp^ '~fp'^p~ and 1]' 'yuj —>■ 'jp'^p". The form factor is proportional to square of 
the residue Cp {cf. Eq.( P^V.6|) ) which is fixed by Eq. (|V.6|) . The results of the calculations 
B{(f) TT^p^ TT^e^e^) = 2.3x 10^^ are in reasonable agreement with the product B'^^^{(j) — >■ 
7r°p°)-B°'''^(p° — ^ e+e") = (1.9 ± 0.2) x 10"^, and similarly for the p+p~ decay mode. 

In case of the p^- and p°-meson decays, the effective vertex has the form 

= -efp^T,erai,v{,d„pldyA^T[^ + d„p1^dyA^Ti~ + d^p^dyA^-K^). (V.7) 



The vertex ( |V.1| ) for V = p*^ and P = 7r° is a part of the vertex ( [V.?] ). Eq.( |V.4| ) remains valid 



for V = p^, p^ and P = tt^ , 7r°. The photon emitted in the decay p tt'-/* is in the isoscalar 
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state, respectively, the transition form factor Fp^T,{t) receives contributions from the uj- and 
0-mesons. The effect of the 0-meson is small due to the OZI rule. The small difference 
between masses of the p- and cu-mesons is beyond the accuracy of the VMD model. Away 
from the p- and tu-mesons poles F^^T,(t) ~ Fp^T^it), whereas F^^.„(t) = F^^T,{t) everywhere. 

B. Decay modes uj r]£^i^ , 7]£^i^ , and (p r]t^l~ 

These decays are treated in the same way as the decays discussed above. The effective 
vertex for the V rj'-y transitions with V = uj, p°, and has the form of Eq. (|V.l| ) where 



one should put P = rj. The decay widths are given by Eq.([VlD with /ip being the rj-meson 



mass. In the simplest version of the VMD model, the transition form factor Fv'^^(t) has 
the form of Eg. ( [V.6|) . For the u- and 0-meson decays, the form factor is determined by the 
p-meson mass, while for the p°-meson decay it is determined by the a;-meson mass. 

C. Decay modes rj' — > cje+e" and rj' — > p^e^e~ 

These decays are of the same nature. The effective vertex for the rj' V'j* transitions 
with V = UJ and p° has the form of Eq. ( |V.1| ). The spectrum of the pair in P ^ Vi~^i~ 
can be obtained from Eq. ([V.4|) with the replacements P ^ V and ^ my. In our case. 



V = UJ and p°, P = rj', the value ^/s = stands for the T^'-mesons mass. The transition 
form factor Fj^i^y{t) has the asymptotics for t oo. It is described within the VMD 
model by the vector meson contributions from isovector and isoscalar channels, respectively. 

D. Decay mode ai(1260) iri+r 

The axial vector meson ai(1260) has quantum numbers I^{J^'-'') = 1^(1++). Its mass is 
above the 0(lO2O)-meson mass, so we do not expect a noticeable effect from the production 
and the decay of such a meson at GSI energies. This resonance becomes, however, important 
at higher energies (see fS^ ). 
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The Oi — > 7r7* transition is an electric dipole El transition. The effective vertex has the 
form 



(5£ai77r = fai-ynid^a"^ - 8 ^a^^) F ^.^'k'' (V.8) 

where F^i, is the electromagnetic tensor. The matrix element for the process ai — > 7r7* can 
be written as follows 

where P is the momentum of the decaying meson, k is the photon momentum, -Fa^^^(t) is 
the transition form factor, -Fai77r(0) = 1. The square of the matrix element summed up 
over photon polarizations and averaged over polarizations of the ai-meson is proportional 
to p*'^{y/s, n, M) + where ^/s = rUa^. The decay width Oi — ^ Tri^i^ takes the form 

rfr(ai ^ 7r£+r) _ p*'(x/i,/x,M) + ^M^ p*{^, fx, M) 
r(ai 777) 0) p*{y/s,n,0) 

X |F,,,.(M2)|^Mr(7* ^£+r)^. (V.IO) 

This expression differs from the analogous expression from Ref. [^], Eq. (6.18), by the extra 
term in the ratio for squares of the matrix elements. Such a term is needed to ensure 
the correct threshold behavior of the s-wave ai iT'-f* decay. 

VI. MESON DECAYS TO TWO MESONS AND £+£- PAIR 

The four-body meson decays proceed either through a two-step mechanism similar to 



that discussed a long time ago by Gell-Mann, Sharp and Wagner [53] in connection to 
the u 7r+7r^7r° decay or through a bremsstrahlung mechanism when a virtual photon 
is emitted form an external meson line. In the first case, the matrix element is a smooth 
function of the photon mass. At small M such that 2me < M, the differential widths behave 
like 1/M, according to Eq.( [ll.l(lD . In the second case, the matrix element is singular at small 
M, as a result of which the differential width increases faster than 1/M. 
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In this Sect., we derive analytical expressions for the double differential widths 

d^T/dsudM of the t] n+TT~i+i-, t]' 7r+7r-£+r,p° ^ 7r+7r-£+r, p± ^ n^n^i+i-, 
fo 7r+7r^£+£~, and — > n^rji'^i' decays, with Si2 being the invariant mass of two 
outgoing mesons. The widths of the other decays are calculated numerically. 

The experimental data are available for the rj Tr+vr^e+e^ decay only, with large errors 



5JI . Nevertheless, the dilepton widths can be predicted using the conventional assumptions 



on a two-step decay mechanism and/or a bremsstrahlung decay mechanism. 



A. Decay modes t] —>■ vr+vr and rj' —>■ vr+vr 

The diagram contributing to the decay r] 7r+7r^7* is shown in Fig. ^. The p^'-fi] vertex 
is given by Eg. (|V. 1| ) . The vertex for the p — tttt decay has the form 

1 a 8 ^ 

= /p..(pj7r+2 n- + p+7r-2 7r° + p^7r°2 7r+) (VI.l) 

where a, (3, and 7 are isotopic indices of the p- and 7r-mesons. The coupling constant 
fp^^ = 6.0 is determined from the p-meson width 

r(p° ^ vr+TT-) = ^jiy{V~s,P,p) (VI.2) 

where ^/s = irip is the p-meson mass. The value fp-j^^^ is in good agreement with the relation 
fpTT-w/gp = 1 which follows from the p°-meson dominance in the pion form factor. In the 



FFGS model the t-dependence of the pion form factor ^^^(t) is attributed to the p- 
meson propagator. It means that the constant fp-j^jr does not depend on the value t when 
the two pions are on the mass shell. 

The radiative decay rj 7r+7r^7 has been measured in Ref. [Q. The 77' 7r+7r^7 decay 
is studied experimentally much better (see p^- pUI and references therein). The last decay 



is more complicated, since two pions occur in the physical region of the p-meson. It was 
shown that an attempt to fit the experimental data within the framework of the cascade 
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model rl — ^ p^^ 7r+7r~7 results in overestimating the events with a vrvr invariant mass 
below the p-meson mass. In Ref . , a nonresonant contribution is added to the amplitude 
and its parameters are fixed by fitting the experimental data. One of the modifications 
considered for the p-meson propagator is the following one 

(si2 - + impYpY^ (si2 - + impYpY^ + Cp (VI.3) 

with a complex value of Cp. 

It was proposed that the existence of a non-resonant term is connected to the chiral box 
anomaly |Q. The anomaly results to a modification of the p-meson propagator similar to 
that of Eg. ( |VI.3| ) with Cp = l/(3m^). The positive sign of the term added guarantees a 
reduction of the number of events on the left hand side of the p-meson peak. The contact 
term describes a production of two mesons in the relative p-wave where a strong resonant 
pion-pion interaction due to existence of the p-meson is present. The final state interaction 
can be taken into account by dividing the bare vertex by the Jost function determined by 
the TT-meson p-wave phase shift. In this way, the p-meson propagator reoccurs, and so we 
get from the contact term the same old p-meson pole term of the earlier models. The net 
effect of the box anomaly is a redefinition of the coupling constant for those models which 
include only the single diagram of Fig. ^ 

We keep the constant term Cp in the form suggested by Ref. |^ and interpret it purely 



phenomenologically as a modification of the p-meson propagator needed to describe the 
distribution of the two-pion events in the rj' 7r+7r^7 decay. The contact term Cp for the 
?7-meson is assumed to be of the same magnitude as for the ?7'-meson. 
The matrix element of the 77 7r+7r^7* decay has the form 

M = -iefp^^fp^nMpelik) {VIA) 

where 

Mp = eraf,u{Pl - P2)r(Pl + P2) aK{ + C p) . (VI.5) 
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Here, pi and p2 are the pion momenta, S12 = {pi + P2Y , k is the photon momentum, 
and rrip and Tp are the p-meson mass and width. Energy-momentum conservation imphes 
P = pi+p2+k. After summation of the squared matrix element over the photon polarizations 
and averaging over the direction of the two pion momenta their in c. m., we get 



-Sp*^{y/s, M)p*2(v^, /i, fl) 



1 



(VI.6) 



S12 - m2 + irUpTp 

with dQi2 = 27Tsin6d6, and 6 being the angle between the momenta pi and k in the c. m. 
frame of two pions. 

The decay width rj —>■ tt^tt^'j* becomes 



We neglect thereby for the p-meson off-shellness in the form factor Fp^ri{t)- The value 
ffryr) = 1-86 GcV^^ (= due to SU {?>) symmetry) is extracted from the p° — > 777 

width with the use of Eq. ([V.3|) . 

The dilepton spectrum in the decay r) 7r+7r^e"'"e~ is described by expression 

dT{r} -K+Ti-tr) = T{ri n+n--f*)MT{-f* ^+r)— — . (VI.8) 

The model for the decay rj' 7r+7r^7* is a direct extension of the model for the 77- 
meson decay. The coupling constant /p-^^' = 1.3 GeV^"*^ (= -^fuj-yn due to SU{3) symmetry) 
can be extracted from the radiative r^'-meson decay width. The two-step decay r]{i]') — >■ 
p^Ti^ 7r+7r^7* is forbidden by C-parity conservation. Constructing the dilepton spectra, 
one should not take into account e.g. channels r]' 7r+7r~7* and rj' —>■ p°7* simultaneously. 
The second one is already contained in the first one. 
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B. Decay mode — > vr+vr 

The diagrams contributing to the decay 7r^7r"7* are depicted in Fig. ^. The first 

two diagrams describe the photon bremsstrahlung. The third contact diagram is added to 
restore the gauge invariance. 

The vertex pyrvr is given by Eg. ( |V1.1] ). The matrix element of the process p'^ -k^ti"^* 
has the form 

M = tef^^F^{M^)er{P)Mr>.e*^{k) (VI.9) 

where 

Mr, = (2Ql - P)r^^^^^ + {2Q2 - P)r^^^f^ " 2^.,. (VI.IO) 

Here, Qi = Pi + k, Q2 = P2 + k, pi and p2 are the 7r+ and n~ momenta, fi is the pion mass, 
and -FttI^) is the pion electromagnetic form factor. The tensor part of the matrix element, 
M.T111 is transverse with respect to the photon momentum, fc, and the p-meson momentum, 
P: 

Mr^k^ = 0, (Vl.ll) 
Mr^Pr = 0. (VI. 12) 

The role of the last equation is the elimination of the coupling of the spurious spin-zero 
component of the p-meson wave function from the physical sector. Gauge invariance requires 
that the form factor entering the third diagram be equal to the pion form factor. 

The square of the matrix element summed up over the photon polarizations, averaged 
over the initial p-meson polarizations and over directions of the pion momenta in the c. m. 
frame of the two pions, has a compact form 

+ ((s - V)(m2 - V + 2si2) + 2si2M2)L(0) (V1.13) 
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where s = = and S12 = {pi + P2)'^- We have also 



and 



F{0 = (VI-14) 
^(0 = ^ln([^) (VI.15) 



C = ^X [^P*{V^2, fi, ( v^, v^, M) , (VI. 16) 

B^ = ^{s + M^~si2). (VI.17) 



The decay width r(p° —>■ tt+tt 7*) is given by 



The integral over the S12 runs from 4/i^to {^/s — M)^. In the case of real photon emission 
(M^ = 0), we reproduce the corresponding expression from Ref. The pion form factor 
can be taken as 

2 

m^p-t- inipTp 

The dilepton spectrum can be obtained with the use of Eq. ( |1I.1CI| ) as 

rfr(p° ^ 7r+7r-£+r) = r(p° ^ 7r+7r-7*)Mr(7* ^ £+r)— — . (VI.20) 

The square of the matrix element at small invariant masses 2me < M and small momenta 
of the virtual photon is proportional to l/u*"^ where u* = B^^j^ is the photon energy in 
the c. m. frame of the decaying meson. The contribution of the region S12 ^ (v^~ 
the integral over the S12 is of the order log^(l/M), and so dB / dM ~ log'^(l/M)/M. This is 
a general feature of the bremsstrahlung mechanism. The square of the matrix element in a 
two-step decay contains no a large parameter. In such a case, dB /dM ~ 1/M. 



24 



C. Decay modes 7r°7r°£+£-, uj 7r°7r°^+r, and uj vr+vr-^+r 

The diagram contributing to the decay p° 7r°7r°7* is shown in Fig. |^. The vertex puj% 
has the form 

= fpuj-n(^Ta,,y{da^rdyp^^-K^ + d^U^d^p'^'K' + d^Urd^p'TT^) . (VI.21) 

The couphng constant /p^^r ~ 16 GeV~^ was determined by Gell-Mann, Sharp and 
Wagner from the u tt+tt^tt^ decay assuming a two-step mechanism a; ^ pvr — »• 



In Refs. |^5| a direct contact term uj vr+vr vr" originating from the chiral anomaly is 



taken into account. The analysis of the u tt+tt^tt^ and pvr, tt+tt^tt'^ decays with 
inclusion of the contact vertex gives a value fp^jn ~ 12 GeV~^. In the contact vertex, one can 
always select a pion pair with quantum numbers of the p-meson. The final-state interaction 
between these pions, as usually, can be taken into account dividing the bare vertex by the Jost 
function determined from the isovector p-wave vrvr-scattering phase shift. In the two-pion 
channels with the p-meson quantum numbers, the ordinary p-meson propagator therefore 
occurs, as a result of which the contact vertex reduces to the ordinary pole diagrams of the 
same old two-step mechanism. 

The coupling constant /p^^Tr can be extracted from the decays p ^ 717 and lj — > 717. The 
one-pole VMD approximation for the transition form factors gives fp^T, = fp^-^gui ~ 12.6 
GeV~^ and /pi^jr = fui-y-nQp ~ 11-7 GeV^^, respectively. In the two-pole VMD model which 
provides the transition form factors with the correct asymptotics these numbers should be 
multiplied by a factor m\/{m\ — m?p). For nix = 1-2 GeV, one gets, respectively, fpi^^, ~ 21 
GeV'i and fp^^ ^ 20 GeV^^. The QCD sum rules |3[ give fp^^ ^ 16 GeV^^ In Ref. g, 
the value /po^Tr ~ 14 GeV""*^ is used. 

We consider contributions to the strong part of the decay amplitudes from the ground- 
state vector mesons only. The two-step mechanism of the u tt+tt^tt'^ decay with a 
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ground-state p-meson in the intermediate state is essentially identical to the mechanism of 
the decays discussed in this subsection. In our calculations, we use fp^,, = 16 GeV~^. 

The vertex u'-fir is defined by Eg. ( |V . 1| ) . The matrix element of the decay p° 7r°7r°7* 
has the form 

M = -iefp^nfu^nF^^^{M^)eriP)Mr^,e;{k) (VI.22) 

with 

■M-Tfi — ^Tcrp\PcrQl\^pK^uQlKku'7y2 2 ~^ ^Tcrp\PcFQ2\^pKpvQ2Kk'u'7y2 2 (VI. 23) 

L^T^ — LJ2 ~ "^oj 

where Qi = Pi + k and Q2 = P2 + k. 
The square of the matrix element 

/r/0 1 P P 

-^MrpM..*-{-gr. + -^){-gpu), (VI.24) 

is averaged numerically over the directions of the photon momentum k in the c. m. frame 
of the pion 1 and the photon. 

The decay widths can be written as 

X f'''-'\tMl^!^imU!^,,^,_, (v,.25) 

J{p+MY 

The integral over the S12 runs from (p + M)^ to (i/i — p)^. The additional factor 2 in 
the denominator occurs due to identity of two pions. The dilepton spectrum in the decay 
p'^ —>■ n^7i^£~^£~ can be constructed with the use of Eg. ( [ll.l(i| ): 

rfr(p° ^ 7r°7r°£+r) = r(p° ^ 7r°7rV)Mr(7* i+r)-^. (VI.26) 

The uj 7r°7r°7* decay is analogous to the decay. The corresponding 

diagrams are shown in Fig. |^. The matrix element has the form of Eqs. ( |V1.22| ) and ( |V1.23| ) 
with interchanged V'jP coupling constants, fuj^.,, ^ fp-y-K, and vector meson masses, ^ 
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The diagrams contributing to the uo 7r"'"7r~7* decay are shown in Fig. ^. The first 
two diagrams (a) and (b) correspond to the common two-step mechanism. The last three 
diagrams occur due to the up mixing and correspond to the bremsstrahlung. The diagram 
(e) restores the gauge invariance of the amplitude. The bremsstrahlung dominates at small 
values of the e~e^ invariant mass. The vertices uop-K and p77r entering into the matrix 
element of the uj 7r+7r^7* decay have the form of Eqs. (|VI.21| ) and (|V.7|) . 



The first two diagrams and the last three diagrams from Fig. ^ are separately gauge 
invariant. These two classes of diagrams contrary to the statement of Ref. |^ interfere: 
Two pions in the two-step part of the amplitude are in an isosinglet state. Due to the 
VMD, the photon in the diagrams (c) and (d) is emitted from the p°-meson. It is equivalent 
to introducing the pion form factor into the tttt'') vertex. The outgoing pion and the p- 
meson with a common vertex form an isotriplet which, in turn, forms a new isotriplet 
with the second outgoing pion. The isospin wave function of the final state looks like 
[7? X [vf X p]] = 7f(7f ■ p) — p{7T ■ vf). The first term corresponds to two neutral pions and can 
be dropped, while the second term describes two charged pions in an isoscalar state. Two 
pions in an isoscalar state from the diagrams (a) and (b) and a photon interfere with two 
pions in an isoscalar state from diagrams (c) and (d) and a photon. The final states are 
identical. 

The interference pattern depends on the production dynamics. The physical cu-meson 
wave function has the form \uj >= \ujq > +e\pQ > where {uq > and |po > are pure isospin 
1 = and 1 states and e is a complex up mixing parameter. The total matrix element is 
given by M = M^P^ + eM^^l The value A^W is the same as in Eq.(VI.22). The value 
is given by Eq.(VI.9). The u 7r+7r^7* decay width represents a sum of three terms 

r{u 7r+7r-7*) = rW(cu ^ 7r+7r-7*) + r^^\u 7i+n--f*) + T^p^^u 7i+7T--f*) (VI.27) 

originating from the vector- meson pole diagrams [p], the bremsstrahlung diagrams [b], and 
from their interference [pb]. It was pointed out in Ref. p2[ that T^^'^^u 7r+7r^7*) = 



2r{u 7r°7r°7*). The bremsstrahlung contribution to the u vt+tt 7* decay equals 
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rM(cj ^ 7r+7r-7*) = |eprM(pO ^ t^+t^-^)^ with |e| = {T{uj 7r+7r-)/r(pO ^ 7r+7r-))V2 = 
3.4 X 10^^ being the p°-meson admixture, as determined from the uj tt^tt^ decay. 
The interference term 

^ 2 1 + ^){-g,.) (VI.28) 

can be found to be 

x „ „ + Cp^(Cp) + (1 - ^(0)/e + (1 - ^(Cp))/g (VI.29) 

where Bp = B^^ + fi^ — and Cp = B^^^/Bp. The value 5^ is defined by Eg. ( |Vl.l7D . The 
interference contribution to the width is given by 



The couphngs /p^^Tr and /p,r7r are determined up a sign. However, the relations f^^-,, ~ fpLjn/gp 
and fpn-n/gp ~ 1 allow to fix the sign of the interference term. The phase 99 of the mixing 
parameter e is a sum of phases ^prod and f decay originating, respectively, from production 
and TTTT decay of the w-meson. The last one is close to 7r/2 (see e.g. Ref. ||6^). We set the 
value ipprod equal to zero. Since the product F.„{M^)* F^^.„{M^) is approximately real below 
1 GeV, the interference can then be neglected. 

The dilepton spectra of the uj — > n^TT^i^i^ decays can be obtained as usually from 
Eq.(II.lO). 

D. Decay modes p r/vre^e^and lu r]TT^e'^e~ 

These decays also proceed through the two-step mechanism. The diagrams are shown in 
Figs. |10] and |Tl]. The matrix element of the decay p rjiT'-f* is given by 

M = efp^J^^^er{P)Mrpe*p{k) (VI.31) 
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with 

1 1 

(cJi — c^2 ~ ^p 

where Qi = pi + k, Q2 = P2 + k, pi is the pion momentum, andp2 is the r^-meson momentum. 
We take the SU (3) symmetry relation fpujirfuj-yri = '^fpp'qfp-i-w into account. 
The average of the squared matrix element, 

/^n 1 P P 

-^MrpM^.*-i-gr. + ^)i-9pu), (VI.33) 

is calculated numerically. The integration over the directions of the momenta of the two 
mesons is performed in the c. m. frame of vr and f]. The decay widths can be written as 



X /''^-"''^P-(v^.v^-^^)P-(v^.^.^-),,., (v,.34) 



where n' is the 77- meson mass. In the form factor F^^^(t), the shift of the w-meson from the 
mass shell is neglected. The dilepton spectrum in the decays p rjT:e^e~ can be constructed 
with the use of Eq. (|irTO| ). 

The uj rjiT^'y* decay is depicted in Fig. |Tl|. The matrix element has the form of 
Eq.( [VI.3lD , its tensor part once the SU{3) symmetry relation fpujnfp^ri = '^fujajrifuj-yn is taken 



into account has the form of Eq. (|VI.32|) , and the 00 — > rjir^'-f* decay width is given by 
Eq.( [VI.33D , with the apparent substitutions fp^ri fuj-yr] and m^^ ^ nip. The dilepton 



spectrum can then be constructed, as usually, using Eq. ( [11.10]) . 



E. Decay mode p+ ir'^Tr^l^i 

This decay is depicted in Fig. The first three diagrams describe the dilepton 

bremsstrahlung. The last diagram corresponds to the two-step mechanism of the dilep- 
ton emission. These two sets of diagrams, each of them is the gauge invariant, interfere and 
we calculate the interference without further approximations. 
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The photon emitted from the p^-meson hne interacts with the electromagnetic current 

jf^{p',P) = + P)^^9ra + KlrQat, ' ga^r/.)] Crb) (VI.35) 

where q = p' — p, ^t{p) is the p-meson polarization vector. The value of A is connected to 
the p"*" magnetic moment. Let's fix first the value of A. 

The magnetic moment of a particle can be calculated from expression 

p = y"rfxi[xx j]. (V1.36) 

The wave function of a vector particle in its rest frame where p = (m, 0) has the form 
er{p) = (0,e). Substituting this expression into Eq.( |V1.35D and Eq.( |V1.35D in turn into 
Eq. dVOeD , we obtain 



p = --^e*se-y. (VI.37) 
Irrir 



where s is the spin operator acting on the space-like part of the p-meson wave functions 
according to the rule SqC/j = —iea/3'ye.y (a,/?, 7 = 1,2,3). In the non-relativistic quark 
model, the p"*" magnetic moment is the difference of the magnetic moments of the u- and 
d-quarks: Pp+ = fiu — ^^d- The proton and neutron magnetic moments are given by the well 
known expressions 

^^r> = (VI.38) 
/X. = (VI.39) 

and thus one can write 

^ 3 

fJ'u- fJ'd= -^{fJ'p- fJ'n)- (VI.40) 



2mp 5 

One gets A = 2.3. The additive quark model predictions are valid with an accuracy of 
10 — 30%. The one-gluon and one-pion exchange currents and the anomalous quark moments 
connected to the violation of the OZI rule, are shown to be important in order to explain 
the deviations of the naive SU{3) predictions from the experimental data ||66|j67[l . 
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The diagram (a) in Fig. 12 contains the pion form factor ^^^(t). A part of the second 
diagram (b) connected to the p-meson convection current contains a p-meson form factor 
Fip{t) which decreases as ~ at t ^ oo. The second part of the same diagram connected 
to the spin current contains a different the p-meson form factor F2p{t) which decreases also 
hke ~ at t oo [Q. It is not clear what kind of the propagator should enter into the 
diagram (c). The gauge invariance requires, within the framework of the considered model, 
that the two form factors Fip{t) and -F,r(^) should be equal. This is, however, a shortcoming 
of the model. It can be overcome on the basis of a more complete theory. The second part 
of the diagram (b) is gauge invariant. The form factor F2p{t) is therefore not subjected to 
additional constraints. The dilepton invariant masses in the p+ 7r+7r°£"'"£~ decay are not 
very large, so the distinction between these three form factors does not exceed a 20% effect. 
We set the form factor -Fip(t) and the form factor entering the diagram (c) both equal to the 
pion form factor. The sum of the diagrams (a) - (c) is then gauge invariant. The form factor 
F2p{t) coincides with the transition form factor F^^T^[t) in the framework of the extended 
VMD model. 

The matrix element of the p"*" 7r"'"7r''7* decay is given then by 

M = tefp^^F^{My^{P)Mr^e;{k) (VI.41) 
Its tensor part consists of the 5 pieces: 

5 

Mr, = J2M',p (VI.42) 



i=l 



where 



1 (pi-p2).(2P-fc), 
(P - - ml 



(VI.43) 



Mlp = A -'•^"7"^V-^-"^^^ (VI.44) 



^ {p,-p2 + kU2p^ + k)p 

(pi + ky-fi^ ' ^ ' 

Mtp = -9r,, (VI.46) 



Pa{P2 + k)xept,py{p2 + k)^,K t\n 

'■Til / / I ; \9 9 ' \ ' 

[P2 + ky - 
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with 

fpunrfcv^nF^jT^i^M ) ^ f pun Fui-yiri'^ ) Afi\ 

The couphng constants fpunr, fu-yn, and fp^n are determined earher (see Table 1). We used 
in Eq. (|VI.48|) the relations f^j-yn ~ fpujir/dp and fpmr/gp ~ 1. It is seen that, while signs of 
the couplings /pTrvr and /p^^^ are a matter of convention, the value 7 is positive and thus the 
sign of the interference term is fixed. 

The first two terms, ( |VI.43|) and (|VI.44| ), originate from the diagram (b) of Fig. |12|. The 



third term ( |VI.45|) originates from the diagram (a), the fourth term comes from the diagram 



(c), and the last term originates from the diagram (d). The second and the last terms 
and the sum of the rest are transverse with respect to the photon momentum: Ai'^^k^ = 0, 
Ai^^k^ = 0, X]i=i 3 4 -^T/x^M — 0- The matrix element, J^rp, is not transverse, however, with 
respect to the p-meson momentum, M-rpPr 7^ 0. A complete theory should yield a matrix 
element which is transverse with respect to both, the photon and the p-meson momenta. 
This is not the case in our effective theory. In the following, we work exclusively with the 
transverse part of the matrix element which is responsible for decays of the physical states, 
and ignore the longitudinal part responsible for a decay of a spurious spin-zero component 
of the p-meson. This is achieved by contraction of the square of the matrix element with 
the polarization tensor 

S.^(P) = -gr. + ^ (VI.49) 
s 

of the massive vector particle, as we did e.g. in Eq.(VI.24). The difference is, however, that 

the term PrPa/s in Eq.(VI.24) can be dropped before calculations. In the case considered 

in this subsection the term PrPa/ s cannot be dropped. 

The square of the matrix element averaged over the initial polarizations of the p-meson 

and summed up over the photon polarizations has the form 

5 

7^ = 5^|>^|^ = 5]7^., (vi.so) 
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where 



P P 

Jr 7--r u , 



7^,, = M\^Mil-{-gr. + ^){-9,u) 



(VI.51) 



Let's now consider separately the different interference terms. 

We perform an additional averaging of the terms TZij with i, j — 1 -j- 4 and i = 1 -j- 3 
and J = 5 over directions of the pion momentum pi in the c. m. frame of two pions. The 
remaining terms 7?.45 and TZ^^ are averaged over the directions of the pion momentum p2 
in the c. m. frame of the pion 7r° and the photon. The decay rate will finally be given 
as a sum of the two integrals over the variables Si2 = (pi +^2)^ and S13 = (p2 + k^- The 
straightforward calculations give 



7^n 
7^22 



4 2si2 + 2s-M2 
~9 (Si2 - sf 

4 



P* (v^,/^,/^)(3 + 



Sl2 



P*\V^, V^, M)-3M^ + 



9 (si2 - sy 

^ 's + M^- S12) (s + S12 -M^){s-M^- S12) 



As 



(s-M^- S12) 



12 



^33 = 



4ss 
1 



12 



3^2 



1 



52(3s + Si2-4/i2 



"^44 
^55 



1, 

4 



7^ 



9 (s23 - m2 )2 



SS23P*^(V^, l^)P*^ {^/s23: 



(VI.52) 



7^12 + 7^21 

7ei3 + Tesi 



T^lb + ^51 
^23 + ^32 



8 A s + S12 - 

9 (Si2 - s)2 S12 



3 (Si2 - 



-(s + S12 - M^){s + S12 + - 8/x2)L(0 

^14 + ^41 = 0, 

47 sp*'(v^, /^)p*'(v^, V^,M)1 + - l)L(C) 



(Si2 - s)B^ 



A 



3 (si2 - 



12ssi2 



-((4/i2-si2)M^ 
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+(-8/x^s - 10ssi2 - 4s^2 + 16At^si2)M2 + (V - Si2)(s - S12)(S + 5Si2)) 



^24 + ^42 
^25 + ^52 



0, 

2 A7 



3 (si2 - s)B^ 



^B'^ilC + 1 - HO) + ^^i2i?^(-l + i^(C)) 



+B^(--sp*W^, M) + (s + M^- s^2)p*\V^2, //))(! - i^(C)) 



7^34 + T^43 
^35 + ^53 

7?.45 + 7?.54 



2 1 



7 



1 



3 5,5, e + c 



-^^12 + /x^)(5i2M^ - \{s -M'- s^mm + cm) 



-i.i2eS.^(-l + L(0 + |(-l + L(C))) 

0. 



(VI.53) 



where 8^ = 3^^ + 11^ — m"^ and ( = Bt^^/ B^. 

The decay width p"*" — > 7r"'"7r°7* is given by expression 



r(p+ ^ 7r+7r"7*) 



q; 



2M2 



X 



(v/i-M)2 4 4 



f%AF^{M')\ 



3 5 



4^2 



ds 



12 



+ / /<-55 "•523 

'(/i+M)2 V'^23 



(VI.54) 



The dilepton spectrum in the p+ — > 7r+7r°^+£ decay can be calculated from equation 



dr{p+ 7r+7r°£+£-) = r(p+ ^ 7r+7r%*)Mr(7* ^ 



(VI.55) 



The decay p — > tt Tr^£'^£ is completely equivalent to the p^-meson decay. 

F. Decay mode /o(980) ir+Tr-i+e- 

The dominant decay mode of the /o(980)-meson is /o —>■ tttt. The effective vertex for this 
decay looks like 
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5C = gf,^Jo7r"7i''. (VI.56) 

The /o vr+vr" decay width has the form 

r(/o ^ vr+TT-) = ±-gl^y(^s,fi,fi), (VI.57) 

with ^/s = rufg. The experimental width r*''*(/o — > vrvr) = |r(/o vr+vr") = 3r(/o 7r°7r°) 
hes within the interval 50 — 100 MeV ||39|. The coupling constant can be found to be 
gf,^^ = 1.6 GeV± 30%. 

The diagrams for the /q 7r"'"7r~7* decay are shown on Fig. O. The matrix element 
has the form 

M = -legf^^M^eUk) (VI.58) 



with 



The tensor A^^ is transverse with respect to the photon momentum. The square of the 
matrix element summed up over the photon polarizations and averaged over the directions 
of the pion momentum in the c. m. frame of two pions has the form 



= A [(M^ - V)F(0 + (2.12 - V - M')m] , (V1.60) 
where su = {pi +^2)^- The functions F{^) and L(^) are defined in Eqs. ([V1.14|) and ( [V1.15|) 



in terms of the parameter ^ which, in turn, is defined in Eq. (|V1.16| ). The decay width 
becomes 

X /■'-^-"'^^•(Av^-^ytv^.^^),,,. (v,.61) 



The dilepton spectrum can be obtained with the use of Eq. (|11.10| ). 
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G. Decay mode ao(980) n-qi+l- 



The dominant decay mode of the ao(980)-meson is ao ttt]. The effective vertex for this 
decay looks hke 

6C = ^ao.^a^vr"ry. (VL62) 
The ao tct] decay width has the form 



r(ao TTT]) = ^c/^„^^p*(yi,/i,/i') (VL63) 



where fi and fi' are the pion and the 77-meson masses. The experimental width r(ao —>■ tti]) 
lies within the interval 50 — 100 MeV. The value of the coupling constant (^qottj? determined 
from the width is given in Table 1. 

Since the 77-meson is neutral, photons can only be emitted from the oq- and vr-mesons in 
an isospin J3 = ±1 state. To the lowest order in a, T{aQ Ti^rj'-f) = 0. Two diagrams for 
the rjiT^'y* decays are shown on Fig. 0. The matrix element can be written in the 



form 

M = -lega^^M^sl^k) (VI.64) 

with 

where pi is the pion momentum, ^/s = nia^ is the Oq- meson mass, and = s. The tensor 
part of the matrix element, is transverse with respect to the photon momentum. The 
square of the matrix element summed up over photon polarizations and averaged over the 
directions of the pion momentum in the c. m. frame of the tt- and ?7-mesons has the form 



- — ^— [(. + su -M' + 2(/.2 - /i'^))L(r) + i?;] (V1.66) 
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where S12 = {pi +^2)^, P2 is the 77 momentum, p| = /x'^, and F{^') and L{^') are functions 
of the parameter 

r = ^./^P*(v^,/^,/^')P*(v^,v^,M). (VI.67) 
The value B'^ is given by 

Bl = hs + M'~ s,2) + - - - (VI.68) 

The a,} ^ ti^t]'-)* decay width becomes 



X r-"" ^P-(v^.^.^y(v^./'.^0,,,^. (v,.69) 



The dilepton spectrum can be calculated using Eg. ( [II. 10]) : 



dT{a+ n+rie+e-) = r(a+ ^ 7r+r/7*)Mr(7* e+e-)— — (VI.70) 
The possible values of M lie in the interval 2m < M < y/s — ji — jj! . 



VII. NUMERICAL RESULTS 

The dilepton meson decays have essentially the same physical origin as the radiative 
meson decays. Eq.( |11.3D allows a direct calculation of the radiative widths within the frame- 
work used for the calculation of the dilepton widths. Thus, the radiative decays provide a 
useful test for the present model. 

A. Radiative meson decays 

The results for the radiative decays with two mesons in the final states are shown in 
Table 2. The p° ii^n^^, uj 7r°7r°7, and r]' —>■ 7r^7r^7 branching ratios are in the 
reasonable agreement with experiment. The branching of the 77 7r+7r~7 decay is higher 
than the experimental value. The 77-meson width is known with an accuracy of 10%. This 
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uncertainty is also present in the calculated branching. The coupling constant /^^^ to which 
this branching is proportional is known with an accuracy of 20%. It is, however, not possible 
to use these uncertainties in order to decrease the predicted branching, since the branching 
for the 1] 7i^7i~e~^e~ decay is lower than the experimental value (see Table 3). The 
measurement of the rj 7i^7i^e^e~ reaction has been done quite a long time ago and with 
rather pure statistics (see Ref. It would be desirable to remeasure this process. The 

experimental branching for the f]' 7r+7r~7 decay is taken from Ref. [^, others are from 
|39| . The results for the p° 7r'^7r°7, p vr°r77, u TT^n^'-f, and u — > 7r°?77 decays are in 



agreement with the calculations by Bramon, Grau and Pancheri where the same model 
for these decays is used. The branching ratios for decays including photon bremsstrahlung 
are calculated for photon energies above 50 MeV. The data on the radiative decays with one 
meson in the final state were used as input to fix the V'-yP couplings. These decays are not 
listed in Table 2. 

The results for the dilepton branching ratios are summarized in Table 3. 



B. p-meson decays 

The branching ratio B{p^ ti^ti' -7) = (9.9±1.6)xl0-^ is almost one order of magnitude 
greater than the branching ratio -B(p*' — > 7r'^7) = (7.9 ±2.0) x 10^^. Despite the fact that the 
first one depends on the experimental cut in the photon energy, one can expect that the decay 
mode p° 7r+7r^e+e~ is important. Indeed, it is seen from Table 3 that the p'^ —>■ 7r+7r^e+e~ 
branching ratio is more than one order of magnitude greater than the p —>■ vre^e^ branching 
ratio and 4 times greater than branching ratio of the direct p° e^e~ decay. The dominant 
contribution comes, however, from the region of small invariant masses of the e^e~ pair, 
which is not seen experimentally. The same is true for the 7r^7r°e"'"e~ decays. For 

invariant masses M > 100 MeV, the number of the e+e~ pairs due to the p^ 7r+7r^e+e~ 
and n'^n^e'^e^ decays increases by 30% as compared to the direct mode, if one assumes 

that the p°- and p^-mesons are produced with equal probabilities. The differential branching 
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ratios for the p-meson decays are shown for the e^e~ channels in Fig. |T5| (a) and for the 
fj^~^fi~ channels in Fig. |1^ (b). At small invariant masses such that 2me < M, the slope 
equals log log M ^ —1 where B' = dB/dM, when there is no bremsstrahlung. Such 
a behavior appears because of the dM/M dependence of the differential branching ratio 
according to Eq.( [ll.l(]| ). The bremsstrahlung makes the M — M'7* widths more singular 
(~ log'^(l/M)) and thus the slope of the other curves is larger. 

In order to get a comparison with the direct channel, we plotted also a weighted — >■ 
e~^e~ distribution according to equation 



dB 1 2MmpTp{M) 



;B{p^i+r). (VII.l) 



dM TT (M2 - m2)2 + mlTp{My 
The p-meson width as a function of the value M is assumed to have to the form 

p*^{M, ^/2, ^/2) ml + p*\mp, ^/2, ^/2) 



V (M) =V ^ ^ ' V ^ P V V u/ , V u/ / , s 



where T p is the total experimental width and ^/sQ = 2p is the two-pion threshold energy. In 
the narrow- width limit, the integral from the left side of Eq. ( |V11.T| ) coincides with i?(p — >■ 



C. Lij-meson decays 

The differential branching for the cu-meson decays are shown in Figs. |16| (a,b). The uj 
decays are dominated through the uj 7T^e^e~ decay mode. The other modes give only a 
small correction to the background. It is interesting that strength of the channel u vr'^e^e" 
is one order of magnitude greater than strength of the direct channel u e^e^ . The e^e~ 
pairs from these two decays have quite different invariant masses. The background appears 
below 500 MeV. Like in the case of the p-meson, we plot also the weighted uj 
distributions according to Eq. ( |V11.I| ) with substitutions nip ^ and Tp F;^ for the 
threshold energy ^/sq = 3p. 
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D. (?!)-meson decays 



The differential branching ratios for the 0-meson decays are shown in Figs. (a,b). 
We do not consider two-meson finals states. The dominant ones could be — > TCTT'-f* and 
(j) — >■ 7r°?77*. The first decay with a real photon in the vr^vr^ mode was observed experimentally 
It has branching 5^^p(0 7r°7r°7) = (1.14 ± 0.22) x 10"^ This decay dominates 
through the /o7 mechanism. One can expect that branching B{(j) -n-^n^e^e^) is two 
orders of magnitude smaller than branching of the radiative decay, and so smaller than all 
0-meson branching ratios listed in Table 3. The — Trrj'-f decay was also measured in 
the SND experiment [HH to give S^^p(0 -> 7T°r]-f) = (1.3 ± 0.5) x 10"^. The branching 
B{(f) TTrje^e^) is also expected to be small. The measured mr'-f and — ^ n'^ri'-f widths 



are significantly greater than predictions of Ref. ||6^ where formation of the intermediate 
/o7 and ao7 states is not considered. At present there is, however, yet no experimental 
evidence for an ao-meson structure in the TT^rj'y decay, and thus it is not clear how to 
interpret the excess in the TT^rj'y mode. Consequently, the calculation of the branching 
ratios for the dilepton TT^rji^i^ modes is still unclear. 

We do not consider decay modes with kaons in the final states. The invariant mass of 
the e"^e~-pairs in such decays does not exceed 50 MeV, so they occur in the region excluded 
by the experimental cuts M > 50 MeV. In Figs. O (a,b), the £^£^ distributions are 
plotted according to Eq.( |VII.l| ) modified for the 0-meson. The dominant two-kaon decay 
channel with y/so = 2 fix where is the kaon mass and also the three-pion decay channel 
with ^/so = Sfj, are taken into account in the M-dependent width T^{M). 



E. r]- and r^'-mesons decays 

The differential branching ratios for the t]- and ?7'-mesons are shown in Figs. |T8| (a,b). 
The 77-meson dominates through the 77 modes. The dominant e~^e~ modes for the 

V-meson are V 'je+e" at M > 250 MeV and r]' n+n-e+e- at M < 250 MeV. The 
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r( ^t^t~ modes display clear structures connected to the uj- and p-meson contributions 
to the ?7'77 transition form factors. The r/ n^n^l'^t' decay is forbidden by C-parity 
conservation, similarly for the r/'-meson. 



F. vr'^-, /o", and oo-mesons decays 



The differential branching ratios for the 7r°-, /q-, and Oo-mesons are shown in Figs. ^ 
(a,b). At M < 500 MeV, the /o-meson decays dominate through the /o n^n^i^i^ mode. 
The four-body ao-meson decay also becomes dominant with decreasing the values of M. The 
dilepton spectra of the /q- and ao-mesons display clear structures connected to the uj- and 
p-meson contributions to the transition form factors. 



VIII. CONCLUSION 

The present work is an attempt to approach the DLS puzzle through including new meson 
decay channels. We studied decay modes of the light mesons below the 0(lO2O)-meson into 
e~^e~ and fi~^fi~ pairs. Besides the direct decays and some Dahtz decays, which are well 
known in the literature, we presented a systematic investigation of the decay modes which 
contribute to the dilepton background and have not yet been taken into account in previous 
studies. These are special Dalitz decays as, e.g. rj' , /q — * •yi'^i' and most of the 

decays to four-body final states. Although many of these processes are found to give small 
contributions to the total dilepton branching ratios, this is not the case for all of them: 

We found that in the p°-meson decays, the dominant contribution to the background 
below 350 MeV comes from the p° Tr+vr^e+e^ decay. The decay modes vr+Tr^e+e" 
and —>■ 7r^7r°e+e~ increase the number of the e"*"e~ yield with invariant mass M > 100 
MeV by about 30% compared to the direct mode p° e^e~ at M < 1 GeV, if one assumes 
that production cross sections for the p^- and p°-mesons are equal. In case of the /o-meson, 
the /o TT+TT^e+e^ decay is in the interval 100 — 400 MeV one to three orders of magnitude 
more important than the next to the dominant /o 7e^e~ mode. In most other cases, 
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the four-body decays give only small contributions and thus there neglection appears to be 
justified. 

The calculated branching ratios can further be used for experimental searches for the 
dilepton meson decays. 

The production rates of mesons in p + p, p + A, and in particular in heavy-ion {A + A) 
collisions differ by orders of magnitude and depend sensitively on the available energy of 
the system. The relative strength of specific channels for different mesons can hardly be 
estimated without complete transport calculations. However, for every meson we can decide 
which channels are more important. The relative strength of e.g. p° —>■ n^i^i" and — *• 
7i^7i~i~^i~ decay modes can be found without knowing the cross section for the p°-meson 
production. 

Since theoretical models for the dileption production have to be adjusted to p + p and 
possibly TT+p reactions before applied to A + A reactions, a detailed knowledge of all relevant 
decay channels and background processes is indispensable. In particular, the lack in the 
understanding of the excess of dilepton pairs below the p-meson peak in the BEVALAC as 
well as in the CERN data shows that scenarios with a reduced p-meson mass can be treated 
at present as a plausible hypothesis only. To draw definite physical conclusions, one needs 
to eliminate possible trivial explanations like those connected to existence of the nondirect 
dilepton decays of light unfiavored mesons. 



After submitting this paper, we got to know on papers by Koch and Lichard ||69 
where a quite complete set of the Dalitz decays (without scalar mesons) and also a few 
four-body channels were investigated. In Ref. |^8[, four-body channels are discussed within 



an approximation suggested by Jarlskog and Pilkuhn |^, which consists in neglecting the 
M-dependence of the matrix elements for emission of the virtual photons with mass M. 
We performed straightforward calculations without approximations and included additional 
Dalitz and four-body channels. In Ref. the decay mode p° n^-ri^e'^e^ is analyzed 
in details. Our expressions (VI. 13) and (VI. 18) coincide with the corresponding formulae 
from the Lichard paper where, however, the pion form factor was neglected. The pion form 
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factor has a small impact on the total e^e decay rate. The dimuon mode, where higher 
invariant masses are involved, is enhanced by about 40%. 
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TABLES 





















ffVl' 


6.03 


1.6 


2.4 


16 


2.33 


0.74 


0.60 


1.86 


1.29 


0.45 



TABLE I. Meson-meson and meson-photon coupling constants. The coupUng constant /ptttt is 
dimensionless, g-ZTrTr and garjw are in units of GeV and fp^T^ ■■■ ffnri' are in units of 1/GeV. The 
meson-photon couphng constants are determined from the radiative meson decays. The extracted 
values are in good agreement with predictions from the SU{3) symmetry. The first three coupling 
constants are determined from the dominant meson decays. 
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Decay mode 



^exp 



7r^7r°7 


4.0 


X 


10-3 




p° — >■ 7r+7r"~7 


1.2 


X 


10-2 


(0.99 ±0.16) X 10 




1.2 


X 


10-^ 




p° — > 7r'^r77 


3.8 


X 


10-10 




uj — 7r'*"7r~7 


3.2 


X 


10-4 


< 3.6 X 10-3 


a; — 7r'^7r°7 


3.1 


X 


10-^ 


(7.2 ± 2.5) X 10-5 


a; — > 7r'^r77 


2.1 


X 


10-7 




?7 — 7r"'"7r-7 


6.9 


X 


10-2 


(4.78 ±0.12) X 10 


77' — 7r"'"7r-7 


2.5 


X 


10-1 


(2.8 ± 0.4) X 10-1 


/o ^ 7r+7r-7 


1.1 


X 


10-2 




^0 ^ ^"^^77 


2.4 


X 


10-3 





TABLE II. Branching ratios of radiative meson decays. For channels with the photon 
bremsstrahlung, the branching ratios are calculated for photon energies above 50 MeV. The ex- 
perimental branching of the rj' — > 'k'^'k~^ decay is taken from Ref. [52], other data are from Ref. 
[35]. 
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Decay mode 


e+e 




^exp 

e+e 






^exp 




input 




(4.48 ± 0.22) X 10"^ 


4. 5 X 10- 


-5 


(4.60 ± 0.28) X 10-5 


p 7r£+^- 


4.1 X 10" 


-6 




4.6 X 10- 


-7 






2.7 X 10- 


-6 




7.0 X 10" 


-11 






5.4 X 10" 


-5 




1.8 X 10- 


-7 






1.7 X 10" 


-4 




6.7 X 10- 


-7 




pO ^ 7r07rO£+r 


7.5 X 10- 


-8 




2.4 X 10" 


-9 




p — >■ 7rr]i~^i~ 


1.9 X 10" 


-12 










oj i+e- 


input 




(7.15 ±0.19) X 10-5 


7.1 X 10- 


-5 


< 1.8 X 10-^ 




7.9 X 10- 


-4 


(5.9 ± 1.9) X 10-^ 


9.2 X 10- 


-5 


(9.6 ± 2.3) X 10-5 




6.0 X 10- 


-6 




1.8 X 10- 


-9 






3.9 X 10- 


-6 




2.9 X 10- 


-8 






2.0 X 10- 


-7 




7.4 X 10- 


-9 






8.7 X 10- 


-10 










(j) i+r 


input 




(3.00 ± 0.06) X lO"'' 


3.0 X 10- 


-4 


(2.48 ± 0.34) X lO-'' 




1.6 X 10- 


-5 


< 1.2 X 10-^ 


4.8 X 10" 


-6 




4) r]e+i- 


1.1 X 10- 


-4 


(1.3l[J|)xlO-4 


6.8 X 10- 


-6 




ri j£+£- 


6.5 X 10- 


-3 


(4.9 ± 1.1) X 10"^ 


3.0 X 10- 


-4 


(3.1 ± 0.4) X 10-^ 


rj — ^ 'K'^'K'~£'^i~ 


3.6 X 10- 


-4 


(1.3 lJ:^)x 10-3 


1.2 X 10- 


-8 




r( j£+£- 


4.2 X 10- 


-4 




8.1 X 10- 


-5 


(1.04 ±0.26) X lO-'' 


r]' u£+£- 


2.0 X 10- 


-4 










ri' 7r+7r-£+r 


1.8 X 10- 


-3 




2.0 X 10" 


-5 




/o ^ 


2.2 X 10- 


-7 




2.8 X 10- 


-8 






1.4 X 10- 


-4 




4.1 X 10- 


-7 




flO ^ ^^+^- 


6.0 X 10- 


-8 




7.4 X 10- 


-9 
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oo ^ 7rr/£+r 4.0 x IQ-^ 1.4 x lO"'^ 

7r° 7^+£- 1.18 X 10-2 (1.198 ± 0.032) x 10"^ 

TABLE III. The integral branchings ratios of the unflavored meson decays to electron-positron 
and muon-antimuon pairs. The experimental data are from Ref. [35]. 
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FIGURE CAPTIONS 



big. 


i: 


Direct decays of vector mesons into electron-positron and muon-antimuon pairs. 


tig. 


o. 

z: 


Decays of scalar mesons into a photon and a virtual photon. 


tig. 


Q. 

o: 


Decays of scalar mesons into a photon and a virtual photon. 


tig. 


4: 


Magnetic dipole transitions V — > P7*. 


Fig. 


5: 


ry-meson decays through a two-step mechanism. 


Fig. 


6: 


p°-meson decays into a virtual photon and two charged pions. 


Fig. 


7: 


p°-meson decays into a virtual photon and two neutral pions. 


Fig. 


8: 


a;-meson decays into a virtual photon and two neutral pions. 


Fig. 


9 


: a;-meson decays into a virtual photon and two charged pions. The first two 



diagrams (a) and (b) are gauge invariant. They correspond to the two-step mechanism of 
the a;-meson decays. The next two diagrams (c) and (d) occurring due to the cup mixing cor- 
respond to the photon bremsstrahlung. . The last diagram (e) restores the gauge invariance 



of the previous two diagrams. 

Fig. 10: p -meson decays into a virtual photon and rj- and tt- mesons. 

Fig. 11: uj -meson decays into a virtual photon and rj- and 7r°-mesons. 

Fig. 12: p='=-meson decays into a virtual photon and n^- and 7r°-mesons. The sum of the 
diagrams (a), (b), and (c) is gauge invariant. The diagram (d) is gauge invariant. 

Fig. 13: /o(980) -meson decays into a virtual photon and two charged pions. 

Fig. 14: ao(980) -meson decays into a virtual photon and 77- and 7r-mesons. 

Fig. 15: The differential branching ratios for the p-meson decays into e'^e~ (a) and 
(b) channels as functions of the invariant mass, M, of the dilepton pairs. The solid curves 
(8 and 9) are the total ratios for the and -mesons. 

Fig. 16: The differential branching ratios for the a;-meson decays into e+e~ (a) and 
(b) channels as functions of the invariant mass, M, of the dilepton pairs. The solid curve is 
the total ratio. The background is dominated through the Tre+e" and np&pL' Dalitz decays. 

Fig. 17: The differential branching ratios for the 0-meson decays into e^e' (a) and 
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(b) channels versus the invariant mass, M, of the dilepton pairs. 

Fig. 18: The differential branching ratios for the 77- and 77'- mesons decaying into e+e~ 
(a) and n'^n' (b) channels versus the invariant mass, M, of the dilepton pairs. The solid 
curves give the total ratios. The narrow structure in the rj' — > 76+6^ decay is connected to 
the a;- meson contribution to the rj' — > 77* transition form factor (see Eq.(IV.5)). The same 
holds for the rj' —>■ decay. 

Fig. 19: The differential branching ratios for the 7r°-, /o(980)-, and ao(980)-mesons into 
e^e' (a) and [x^ ijT (b) channels versus the invariant mass, M, of the dilepton pairs. The 
solid curves No. 6 (a) and No. 5 (b) give the total branching ratios, respectively, of the e+e~ 
and ji^ji~ modes of the /o-meson. The structures in the /q- and Oo-mesons Dalitz decays 
are connected to the uj- and p-mesons contributions to the /o and oq transition form 
factors. 
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